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ABSTRACT 



A preliminary computation of the Isgur-Wise universal form factor using a lattice 
formulation of the Heavy Quark Effective Theory (HQET) is described and 

compared with the recent data from ARGUS and CLEO on B -^D * I v decay. 



This letter describes a lattice computation of the Isgur-Wise function, also called the 
heavy quark universal form factor. The calculation uses the lattice implementation of the heavy 
quark effective theory developed recently by the authors ^ The calculation is formulated and 
carried out entirely within the Isgur-Wise limit, in which the masses of heavy quarks are treated 
as infinite. We show the first lattice simulation results based on this formulation, and compare 

them to the recent data from ARGUS and CLEO on B -^D I v decay, the actual process that 

comes closest to realizing the Isgur-Wise limit. Two other lattice calculations of this decay form 
factor have recently been published^'^. They treat the heavy quarks as Wilson fermions with a 
small hopping constant, but do not implement the heavy quark Umit. 

We briefly review the lattice formulation of the heavy quark effective theory. After Wick 
rotation to Euclidean space, the classical velocity, which is only a parameter, is complex when 
expressed in terms of Euclidean four-momenta: 



with the boundary condition that S^'^'ix) vanishes for negative x^. In the presence of gauge 

interactions, the derivative is replaced by a gauge covariant derivative.In order to simply 
incorporate the boundary condition on the lattice, we represent the time derivative in the reduced 



V = {ivQ,v) 



(1) 




The reduced propagator for the heavy quark satisfies the Euclidean space equation 



ivDS^''\x) = -iv[d - igA{x)^S^''\x) = 5 (x) 



(2) 
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Dirac equation by means of an asymmetrical forward difference, but use a symmetrical centered 
difference for the space derivative. Tlie lattice equation for the reduced propagator is: 

3 • (3) 

~IV 

+ [C/^(x,x-HA)5^"\x+A,y) - [/^(x,x-A)5(^V-A,y)] - Hx,y) 

This is solved by simple forward recursion; no iterative techniques are needed. A symmetric 
time derivative could also be used, at an increase in computation. 

Even for non-zero v, only trivial computation is needed to construct the heavy quark 

propagator. This makes it entirely feasible to compute three-point functions, and a careful 
choice of heavy-light particle wave functions will improve results as in other lattice calculations. 

The lattice heavy quark theory has some unfamiliar features which are worth noting. The 
symmetrical first difference results in fermion doubling. However, unlike the usual fermion 
situation, the contributions of the secondary modes automatically vanish in the zero spacing limit, 
without the necessity of adding a Wilson-like term. The reason is that they always occur in 
association with either a light Wilson quark or a gluon whose momentum is at the edge of the 
Brillouin zone, and those modes have energies on the order of the inverse lattice spacing. This 
may be seen in perturbation theory^. Note that to leading order in the inverse heavy quark mass, 
there are no closed heavy quark loops. 

Another feature of the lattice heavy quark theory is that for residual momenta oriented 
oppositely to the classical velocity, the propagator grows with Euclidean time. The reason is that 
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the falloff is governed by the difference between the actual energy of a mode and that of the 

zero-momentum mode, which is 

lim^M^ + (Mv+pf - M^l = (4) 

in the heavy quark limit. When this is negative, the "falloff becomes growth. Conversely, at 
fixed Euclidean time, the heavy quark propagator grows in some directions in momentum space. 
The modes of the light quanta always compensate for this growth in simulations. In perturbation 
theory the contour in the Euclidean energy plane is determined by the boundary condition that 
heavy quarks propagate only forwards in time. The use of light staggered fermions would be 
problematic, since they have poles near the edge of the Brillouin zone. 

The paradigmatic heavy quark process is B — >D * / v , the weak decay of one meson 

containing a heavy quark into another. In the heavy quark limit, there is only one form factor 
that describes B meson decay into both D and D* mesons, and it is a function only of the 
classical velocities''^, 

B Dlv , D*lv 

(5) 

{Dv'\cy^^(l -iy^)b\Bv) ~ ^(vv') 

All the dynamical information is contained in the function ^ , known as the Isgur-Wise function. 

It is a function only of the classical velocities, and is normalized to unity for forward decay 
because it also describes the conserved vector current taken between states of identical heavy 
mesons. 



4 



^(1) = 1 



(6) 



For the process B — >D * / v the structure of the matrix element is 



{D*v\e\cy^^(l-iy5)b\Bv) = M^Mj^. C^b^i^'^') r^^(vy,e) 



where F is a kinematic trace of Dirac matrices and projectors, and the normalization constant 
comes from the renormalization point dependent matching of the heavy quark effective field 
theory to QCD with quark masses much larger than the renormalization scale. It has been 
calculated to leading log approximation^: 



C, 



cb 





6/(33 















a(vv') 



(8) 



where A'^ = 4 is the number of "light" flavors below the heavier heavy (B) meson, and 

8[vv'r(vv') - 1] 



a(vv') 



33 - 2N^ 



(9) 



= 3 is the number of "light" flavors below the lighter heavy (D) meson, and 



r(vv') = 



In 



(vv'+ ]l(vvf - 1 



(10) 



(vv'r- 1 
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The renonnalization scale |j dependence of the second factor in Eq. (8) is compensated 
by a scale dependence in ^, in just such a way so that physical amplitudes do not depend on |a. 

The amplitude for B — >D * / v is contained in the reduced three-point function 

G^'«''^\xo,Zo,yo) ~ {Trs{y,x)Tj,S^'''\x,z)TjS^''^\z,y)T^ ^^^^ 

The quark flow in this process is as shown in Figure 1. The residual momenta are eliminated 
by summing over all lattice sites on the Xq and Zq time slices. The Isgur-Wise function may be 
extracted from the three-point correlation by dividing out the normalized single meson 
propagators. 

1/2 1/2 M^b^^d). _ ^ ^ 

q(vv) = lim U^j 



The wave function renormalization constants are computed from the normalization of the 
asymptotic heavy-light meson propagators. Their evaluation, and the attendant numerical errors, 
can be avoided by taking advantage of the normalization of the Isgur-Wise function in the 
forward direction. We use three-point functions with both classical velocities set equal to v and 
to V as normalization, and the expression from which we extract the universal form factor is 



^(vv)|2= lim ^^"^^^(^0>>-0>^o) ^^"^^(^O-^O-^o) (13) 
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This forai has the advantage that the overall normalization of the three-point function drops out 
completely. Only those renormalizations that depend on the classical velocity need be included. 

We evaluated the three-point functions using the lattices made publicly available by 
Bernard and Soni^*^. There were an ensemble of 16 lattices of size 16^ x 24 with lattice coupling 
P = 5.7. The light quarks were represented by Wilson quark propagators with hopping constant 
K = .164. The heavy quark propagators were evaluated for four values of each component of 
the initial and final quark classical velocity: 0, .25, .50, and .75. That is, there were 4^ = 4096 
combinations of and Vq. Of course, the many that were lattice rotations of each other were 
averaged to improve statistics, but even so there were more than a hundred different values of 
Vg-Vj-f, lying in the interval [1,1.822]. 

To improve the overlap of the meson wave function with the ground state, we applied 
multiple single -link smearing steps. The light quark propagators were evaluated with point 
sources, and we applied a "Gaussian" smearing of the form 

(1 + aH)" (14) 
to the heavy quark sources. H s a symmetrical single link displacement operator. 

3 

Hiy,x) = [Uix+i,x)by^^i + Uix-i,x)dy^^_i] (15) 

The parameters used were n = 0, 5, 10 and a = 4.0, following the study by Alexandrou et al. 
for static quarks In the present case, we found that n=5 gave statistically cleaner signals than 
n=0, but that n=10 was not an improvement over n=5. 
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In order to project out the mesonic ground states from the three-point fiinction, one either 
needs excellent wave functions or large Euclidean time separations between the slices on which 
the heavy meson is created, emits a weak current, and is detected. Unfortunately, the statistical 
noise in the simulation overwhelms the signal when those separations are as large as 4 units of 
Euclidean time. The largest separations that provide a signal are 

3^4 - ^4 = 3 

At 2 units of Euclidean time separation, which corresponds to a greater contamination with higher 
mass states than at 3 units, the computed form factor is flatter. This presumably reflects 
Bjorken's sum rule, which states that the sum of the squares of the Isgur-Wise functions to all 
final states is 1^^. 

The results of the simulation together with a comparison to the data on B — >D * / v from 

both ARGUS and CLEO are shown in Figure 2. The agreement with the ARGUS and CLEO 
data with the simulation using 3 units Ar of time separation is obviously rather spectacular. 
However, to know that such a results is reliable we must, at the very least, have available several 
values of Ar and see that the result is stable at the largest values. Since there is a large change 
in the Isgur-Wise function between A? = 2 and A? = 3, and since the signal has degenerated into 
noise by Ar = 4, we simply have no way of knowing whether this agreement is significant or 
fortuitous. 

A simple quadratic fit to the simulation results at A? = 3 gives for the slope of the 
Isgur-Wise function at the origin: 
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^'(1) = -.95 



(17) 



For the reasons just stated above, we regard the numerical results presented here as 
preliminary. There are corrections to apply and improvements to make, beyond the usual resort 
to larger ensembles of larger and finer lattices. In the present calculation we have neglected both 
the matching of the continuum HQET to QCD with quarks whose masses are far above the 
dynamical QCD scale, and the matching of the lattice HQET to the continuum version. One 
estimate of the latter is that it is a small effect^, but it is necessary to study it systematically. 
A significant improvement would be to find operators for the B and D mesons that match the 
states very well. While the "Gaussian" smearing we have used improves the signal somewhat, 
far more systematic techniques have been developed for application to static quarks , and they 
should be extendable to lattice heavy quarks with finite classical velocity as well. A further 
enhancement would be the use of lattices generated with so-called improved actions, and 
propagators from which the leading finite lattice spacing errors have been removed. 
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FIGURE CAPTIONS 



Figure 1. The heavy lines represent b and c quarks, the light one a u or d quark, and the wavy 
line the emission of a weak current. 

Figure 2. The Universal Isgur-Wise Form Factor plotted vs the classical velocity product v- v'. 
Five wave function smearing steps were used in each simulation. The top graph corresponds to 
2 time steps each between the initial heavy meson source and the emission of the weak current 
and between the emission of the current and the final heavy meson sink. The middle graph 
corresponds to 3 time steps each, and the bottom one to 4 time steps each. The error bars on the 
simulation points are statistical only, and include averages over all configurations related by 
lattice symmetry. The solid lines are simply quadratic least squares fits to the simulations to 
guide the eye. The open squares are measurements from CLEO and the open diamonds are from 
ARGUS. 
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